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Abstract 

In Hilbert space, a linear source-to-flux problem in the critical 
(zero eigenvalue) limit is ill-posed, but regularized by a constraint 
on a linear functional, fulfilled by tuning some control variable. For 
any exciting perturbation, I obtain, by spectral decomposition and 
perturbation theory, the regularized flux and the regularizing control 
variable non-linear responses. 

May the exciting perturbation be obtained, inversely, from observ- 
able responses? Yes, in some cases, from the existence of a weight 
scale, a perturbation series, determined by recursion relations, in- 
volving well-posed source problems, and the possibility of obtaining 
this weight scale from observables of both the unconstrained and con- 
strained systems. 
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I Problem setup 

In Ti, a real Hilbert space, with the scalar product (.,.) and the norm 

II ^ 11= {ip,ipy/'^, let L : — i> Q = L(JF) be a linear (additive and ho- 
mogeneous) endomorphism (JF, Q G Ti. are linear manifolds). 

I consider the source problem, finding, for some source Q E Q, the fiux 
ip E J-' such that 

= Lilj + Q. (1) 
I assume that the 'fiux-to-source' operator L is invertible (on the left), hence 

V' = $(L,Q) = -L-iQ. (2) 

It occurs that the fiux operator $, defined in (||), is linear in Q, but 
non-linear in L. This non-linearity is formally eliminated by the change of 
variable L H-i> L^^. However, expressing an inverse operator is often difficult, 
in particular, when the operator L has an eigenvalue close to zero (in this 
case, the non-linearity is somehow traded for a singularity). Therefore, I will 
keep the variable L, and the non-linearity L $. 

The small eigenvalue critical limit of the source problem is conveniently 
defined by the following hypotheses: 

• L,Q are functions of some scalar variable z (varying over an interval); 

• a{z) is a single eigenvalue of L{z) associated with an eigenvector ip{z); 

• is a critical value, such that 

lim (t{z) = 0; (3) 

• for all z 7^ Zc, L{z) has an inverse on the left L~^{z). 

Even though the source problem has no more than one solution for z ^ z^, 
the source problem (|I]) is still ill-posed in the critical limit z — > Zc, because 
$ is infinitely sensitive to perturbations of Q (essentially: z ^ <^{z)~^ has 
no derivative at Zc). 

1 take for regularizing condition the permanence of some linear continuous 
'gauge' functional of fiux. This will be achieved by tuning the variable z, 
thus appearing as a control variable. Equivalently, in the language of system 
analysis, a high-gain amplifier can be hardly used without feedback control, 
in order to maintain a reasonable output level (not too large, not too small). 
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As the Hilbert space is self-dual §9-5, p. 264], the regularizing condition 
expresses as the constraint 



The system parameters are T = (L, Q, Q^), all functions of z. T may be 
excited by some perturbation 6T, actually driven by some exciting variable 
e, independent of z. 

Primary observables are the scalars outputs of the system, here 

• the exciting variable e, 

• the control variable z, 

• the gauge output R. 

Secondary observables are parametric expressions thereof, independent of 
the system parameters T, 6T, for example. In R, or the partial derivative 



Observation is useless, without an interpretation relation, between T, 6T 
and observables. The main issue of the paper is to prove the existence of an 
interpretation relation, of the form 



where Zi, Z2 are scalar weight functionals. The weight Zi is obtained by 
processing the primary observables through Z2 and (|^). I will explicitly 
construct Zi, Z2. 

Ideally, the sources are unexcited {0 = 6Q = SQ'^) and the exciting per- 
turbation is uncontrolled {6L does not depend on z). Solving (H) for 6L, 
if possible, is the principle of a direct operator perturbation measurement 
method. 

2 The source system and the critical Hmit 

2.1 Unconstrained system definition and properties 

The source system is represented by (|1]). As the constraint is not used, the 
control variable z is left constant and implicit ('asleep') until the critical limit 
is considered (section p.3|) . In the language of system analysis, the control 
loop is open. 



R= {Q\<l>) = Ro,Q^ en. 



(4) 



(dR/dz),. 



Zi{T,6T) + Z2(prim. observables) = , 



(5) 



2 THE SOURCE SYSTEM AND THE CRITICAL LIMIT 



5 



The gauge output is 

RiT)^{Q\^L,Q)),Q^ en. (6) 

I assume that L has an adjoint L"!" : JF"!" ^ Qt = L{J^'^) and G . The 
adjoint source system is 

= lV^ + Q^. (7) 

I assume that L^ has an inverse L'^~^ : <^ ^ T\ solving the source 
problem (|^: 

VQ^ G Q\ = gt) = -Lt-iQt = $(2^t^ Qt). (8) 
I introduce a restriction of the scalar product, noted (.|.), to 

(Qt X ^) U (^t X Q). 

Considering, adjointly, the source as the gauge functional of the adjoint 
source system, the adjoint gauge output is 

i?t(Tt) ^ ($t(L,Qt)|Q), Tt ^ (L,gt,g) (9) 
One shows that L^~^ is adjoint to L~^: 

Lt-i = L-i^ (10) 
Using ( [T0| ) in (|]) yields the reciprocity relation 0, 

i?(T) = (gt|$(L,g)) = (<|.t(L,gt)|g) =i?t(Tt). (11) 

From (H, the adjoint flux $"''(L, Q^), considered as a functional of Q, 
produces the gauge output. 

The linearity of the fiux-to-source operator and the gauge functional im- 
ply scaling laws: 

Va G R*,aRo = (aQ^|$(L,Q)) = {a^^{L,Q^)\Q). (12) 

Considering the source and the flux as the physical quantities, while the 
adjoint flux and the gain factor a are mathematical and operational com- 
modities, ([T^) shows that the physical quantities are independent of a, which 
corresponds to what physicists call gauge invariance 0, §VII.5]. 

The unconstrained system is represented by the parametric expression 

f/(T)^(i?,$,<i.t)(r), (13) 
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mapping the independent variables, or parameters, onto the dependent vari- 
ables, the gauge output and the fluxes. With ([T3|), (|T2|) becomes 



Va G M*, U{L, Q, aQ^) = {aR, $, a<^^). (14) 

By definition (Q), and adjointly, the direct and adjoint source systems are 
uncoupled: 

2.2 The spectral decomposition 

I assume that has an eigenvalue cr"'", and that the associated eigenspace 
is of dimension one and eigenvector (yj"'', 7^ 0. Hence, the eigenvectors 
can be and are normalized, so that 

(V'V) = 1- (16) 

(As cr 7^ 0, (y9 G Q, and the use of the restricted product above is correct.) 
One shows (classically) 

= (at-a)(y.t|y,), 

hence, from (|16|), = o . One shows that T = v?^"*" fl ^ is stable for L, and, 
adjointly, jFt = (p-*- n JF"!" is stable for L"!". 

is the hyper-plane, orthogonal to i/?^, a subspace of codimension one. 
From (|16|), ^ y^"''-'-, so that Ti = ® v^^^ and one shows 

which are equivalent to 'closure relations', 
where 

• ifip"^ stands for the tensor product ly^Xv^^l, 

• stands for the tensor product \ip'^){ip\, adjoint to the latter, 

• vf is the projection, in JF, on JF, along ip, 

• tt"'' is the projection, in T'^ , on T'^ , along adjoint to the latter. 
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Condensed notations will be necessary: 

\/ip eJ^,^ = n^, eJ^\^ = L = nL,lJ = n^Ll 

Applying tt to (|l]), 

(} = Lij + Q = L^ + Q, (18) 
which is solved for the harmonic flux: 

$ = -L^^Q, (19) 

where L has been restricted to the stable subspace JF. Applying {^p^ to (|l]), 

(^t|$) = (20) 



a 



Combining (0, l^, |T|), 



$ = -i^v^-L-g, (21) 
a 

R = -i^i^^ + (Qt|_Z-g). (22) 
cr 

On (^TD appears the (unconstrained) source problem critical ill-posedness. 
Eliminating a from (E2|) and inserting the result into (pT]), 



_ (gt|-L-iQ) (gt|j,) 

= ^ - W ^ ^ 

it(l — U) 

u; is the 'harmonicity'. The negative source-to-flux operator is 



Some useful relations are, adjointly, 



^ (-Lt gt|g) ($t|g) 

" = — ^ — =mQ)=^^ 



(27) 
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and 



TT^ = TT, 

7fL = Lit, 



L-i = L~\ 
{Q^\-L-'n = (<l>t|^, (28) 



and adjointly. 



2.3 The constant gauge output critical limit 

I show up the (smooth) scalar control variable z. Acknowledging the func- 
tional character of the parameters T, (|T^) must be corrected to show chained 
functions: 

UoT = {R,^,^^)oT. 
I assume that observables are bounded in the critical limit (H). Consid- 



ering moreover (|2^ , P5| , pTl) , the direct flux or the adjoint flux must diverge. 
I assume 

{ip^\Q)ip{z) = Oia{z)),z^z,, 
3mi,V^,0 < mi <|| {Q^!f)lf\z) ||, (29) 

(while the source is turned off, the gauge is not) so that the direct flux 
does not diverge and the adjoint flux does, though only by a scalar factor. 



Following the comments of section on gauge invariance, the direct flux 
boundedness is physically motivated. 

I assume spectral separation, uniform in z: 

3m, Vz,V(^ G ^, II ^ 11= 1),V(5^ G ?t, II ^ 11= 1), 

< m < min(|| L{z)ij \\, \\ D{z)^^ ||). (30) 

Consequently, any eigenvalue of L has a modulus greater than m. Only 
the eigenvalue a{z) gets close to zero. a{z) and ip{z) are the 'fundamental' 
eigenvalue and eigenvector (and adjointly). The maximum of all m in ( pO] ) is 
the 'spectral gap'. The harmonicity and the ratios of the second terms over 
the first terms in the r. h. s. of ( |25| , p^ ) are 0{a{z)/m). The constant gauge 
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output critical limits are 

ujoT{z) = 0{a{z)/m), (31) 

3 Perturbation theory with constraint 
3.1 The constraint operator 

The constraint consists in tuning the control variable z to the 'balancing' 
value z, such that the gauge output remains equal to a reference value Rq, 
which is a new independent variable, while ^ is a new dependent variable. 
I assume that the problem has a unique solution: 

3\z{T,Ro),R{T{z)) = Ro. (32) 

Remark 1. As T is a function (of z), z{.,Rq) is a (non-linear) functional 
( of the function T ). 

For all function t of z, like T, is defined the constrained value 

t_{T,Ro)=t{z{T,Ro)). 

The (linear) 'constraint operator', 

t \^ t = t o z, (33) 

transforms a function of z into an operator on {T,Ro). (|33D is formally the 
commutation relation 

J = t_. (34) 

For all operator u on T{z) (like Id, U), u{T) = uoT = u{T_) >u. The symbol 
> introduces abbreviations (to be used with care). 

The constrained system is thus represented by {U_,z_), where $ should 
be evaluated from (|25|), where R — > Rq. On the constrained system, (|T^) 
becomes 



Va e M*, iU,z)iiL,Q,aQ^),aRo) = ((«^, «$^), £)((!>, Q, Q^), 

(35) 
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The adjoint source although it does not appear in (|I]), nevertheless 
affects the flux, because the flux-to-source operator depends on z, that de- 
pends on Q^] the constraint couples the direct and adjoint source systems; 
as opposed to ([T5|), 

7^ 0, ^ ^ 0. (36) 



agt ^ ' dQ 

3.2 Functional variational definitions 

The perturbation of any vector x is noted, as usual, 

{x, Sx)* = X + 6x > X* . 

X, X* are respectively the reference and perturbed values. The deflnition is 
extended to any parametric expression f{x): 

ifix),6f,6xy ^ if + 6f)ix + 6x)>nx*)>fixy, (37) 
Sifix),5f,5x) ^ if + 5f)ix + Sx)-f{x)>Sifix))>Af. 

In particular: 



Af = f{x*)-f{x) + 6f{x*), (38 



e. g., for / : X I— ^ y.x, where '.' stands for a bilinear product, 

S{y.x) = y.6x + 6y.x*. (39) 

I assume that the source system remains of the form (0, hence the 
'law invariance' statement (on the operator U, not its value!): 



6U = 0[ (40) 



and, consequently, S{U_,z) = 0. On the contrary, (the function) T may vary 
by (the smooth function) ST. Control is linear if and only if the second 
derivatives vanish: 

= T" = 5T". (41) 
Deflnition: control is remote (from excitation) if and only if 

6T' = 0. (42) 
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A perturbation of the unconstrained system U(T{z)) is, from (0, |iO| ), 

U{T{z))* = U{T{z) + AT) > U\ 

T{z*) - T{z) + 5T{z*) > AT. (43) 

Taking advantage of (^5l), the perturbed constrained system U_{T*,Rq) is 
rescaled, according to the gauge transform 

(Qt*,i?*,$t*) ^ a(Qt*,i?*,$t*)^^ = (44) 

Kg 

SO that 

Ro = Ro- (45) 

Whenever the gauge output is constrained, it is taken as constant, without 
loss of generahty, which allows to drop the independent variable Rq. 
From (pHD, for all operator u on T{z), 

\/{u,6u = 0), Au = u{T*) — u{T) > u* — u> 6u, 

_A = (46) 

the result of constraint and perturbation on u(T{z)) does not depend on the 
order of these operations. In particular, with u G {ld,U}, considering (^OD, 
AT = 6T, AU = 6U_. I will also often use 5T, and one must be careful that 

Remark 2. 6T^ST = AT. 

3.3 Perturbative expressions of the gauge output and 
the fluxes 

Using (l39| ) on (|l|, H), and noticing that U depends on T, that depends on z, 

AR = (Agt|<l>*) + (gt|A$), (47) 
-LA<I> = AL$* + Ag. (48) 

Muhiplying (H) by ($^| and using 

(Q^|A<I>) = (<I>^|AL$* + AQ). (49) 

Combining (P, |9]), 



AR = dRriAT) + ($^|ALA<I>) + (AQ^|A$), (50) 
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where dR is the differential of R{T): 

dRridT) = {<^^\dL<t> + dQ) + {dQ^<t>). (51) 



The harmonic part of A$ is obtained by projecting (^Sj) on the harmonic 
subspace and inverting: 

A^= -L-\AL^ + AQ)-L-^ALA^. (52) 

I abbreviate (pc><^ (and adjointly). The fundamental amplitude (v9^|A$) is 
determined, independently of the fundamental eigenvalue, from (|47| , p!7|) : 

AR = (gt|A$) + (AQt|$) + (Agt|A$), 
(Qt|A$) = (gt|y,y,t|A$) + (gt|A$). 

From (|52| , ^Tj) and after rearranging, 



AR = dRriAT) + {Q^<f<f^\A<l>) + {<^^ALA<!>) + {AQ^\A<^), 
f = {L,Q,Q^). (53) 



Remark 3. T ^ (Z,Q,gt). 

The constraint is applied to (|^, |52|, ^3|). In particular, from (0, ^5]) 



= AR = 5R. (54) 

= (ii?T(5T) + + ((5Qt|5$^^ (55) 

= (5L $ + (5g + (56) 
-(Qt|^^|5$) = (ii^Tl-^f) + + (5Q1-|5$). (57) 

In agreement with (0), not only the direct source, as expected from (P, 
but also the adjoint source, affects the flux, by ways of (^S]). The similarity 
between the r. h. s. of (|55| , |57D helps calculations. 

( |55| , p6| , [57| ) do not yield perturbed quantities in a closed form; but, with 
ST = C(e),e ^ 0, they yield 5$ up to C(e^): they are 'perturbative' 0, 
which allows to solve the perturbation problem by perturbation theory ^ , 
as follows. 
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3.4 Perturbation series of the fluxes and the control 
variable 

Perturbation theory transforms f{x)* into a function of a perturbation vari- 
able e, according to 

Sx eSx. (58) 

This is apphed to the constrained perturbed system. The response is sought 
as a power series of e: 

oo 

{U,z){T + e5T)t>{U,zy = J2{U,zUT,5T)e^. (59) 

n=0 

By definition, {U_,z)n{T, .) is homogeneous: 

Va e M*, {U,z)n{T,a6T) = a"{U,z)n{T,6T). 

I do not study the convergence of the perturbation series: I only find 'formal' 
power series solutions. 

Applying (H, g,||) on T{z), 

T* = T{z*) +eST{z*). (60) 

I assume that (T,6T){z) have power series at z(T). ( |106| , |107| ) are used 
to develop in powers of e the functions of z* in the r. h. s. of (|60D 



T* = + ^n-i^ + T,iz„ . ..z^_,) +6T^_,iz„ . ..z^_,))e\ (61) 



n=0 



where the term n = is correct from ( p.09| ) and a variable shift on T, so that 
= ^0 = = 5T__i- 

The power series (^, ^1]) are introduced into (^, 15^). By convention, 
^0 = 0, and I define, for all endomorphism D on T, and also for D = S, 

Vp>0,(DT)op = {DlM^±DQ^\^^), 

{DT),o = {^,\DL^±DQ), 
{DT) 

Op+pO = 

{DT) = {DT)oo = dRTiDT), (62) 

Vp,g>l,(Dr)p, = {^p\DL^)- (63) 
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By identification at order e": 

Vn > 1, -zjT') = z^^^T') + (T„ + 5T„.i) + 

XI Un-p(^')op + ^„_p-i('5T')op + (r„_p + (5T„„p_i)op), (64) 

l<V<n-l 



U^-pH + z^-p-iMI + L^-p + SLn-p-A, (65) 



l<p<n-l 



Vn > 1, -(gt|^t|$j = ^^(r') + z^_,{6f') + (T„ + 5T„„i) + 

U^-p{T')op + z^^p-i{ST')op + {fn^p + 6fn^p^i)op). (66) 

l<p<n— 1 



Furthermore, using (0), ^„ is eliminated from ( p3| , |5B| ) to obtain recursion 
relations. After some rearrangements: 



Vn > 1, -L i„ = ^„_,(-i^(L'$ + Q') + + 5g') - 

^^^^i^^(r$ + g') + (4 + 5L„„i)$ + £^ + 6Q^_^ + 
E (^„-,(-^(r$+g')+l.'$,)+^„-,-i(-^(r$+Q'^ 



l<p<n-l 



{T') 

{Tn-p + STn-p-l)op , -/ ^ ~l 



(T) 



(L $ + Q ) + (L,_^ + 5L„_^_i)$,), (67) 



^._i((5T') - (5T')|^) + ((T„ + 5T„_i) - (T„ + 5T„-i)|^) + 
J] Un-p((^')op - {T')op^) + z^^p-i{{ST')op - {6T')op^) + 



l<p<n-l 



{{Tn-p + STn-p^l)op — {Tn-p + ^Tn-p-l) Qpj^)) . (68) 
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Finally, the coefficients {U_, z)n can be constructed, up to any order, from 
the reference flux-to-source operator, source and constraint, and their ex- 
citing perturbations, with the recursion relations ( p^ , |57| and adjoint 
relations. First and second order results are 

(70) 



-z,{r) = z,i{5T') + (T%0 + (T,) + {5TU ■ (71) 



In ( [7lD appear the various contributions to control variable response 
non-linearity: (T2), {ST') represent, respectively, the intrinsic control non- 
linearity and the cross effect of control and excitation, vanishing if control is, 
respectively, linear (^) or remote (^). (T')oi, (5T)oi represent control and 
exciting perturbation shielding, as they involve the flux shielding coefficient 
$2; ^ global consequence of both exciting and control perturbations. 

Remark 4. Because of shielding, the control variable response (e ^ z_) may 
he non-linear, notwithstanding linear and remote control. 



i = z,{-^^l^p^cil^ + g') + + + 5^ 



- (QW^I^2) = zMSf') - (5T')|^) + {{f'U - (T')oi|^)) + 

((^2) - (T2)|^) + - (5T)oi|^). 

For linear control, using ( |111| ), (^^, ^7\, |BS|) become 

Vn > 1, -z^{T') = z^_,{6r) + 6nAST) + {6T)o,n-i + 

J2 Un-p{T'K + Zn-p-l{ST')o,), (72) 
l<p<n-l 
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Vn > 1, -L = z^_^[J^CL^ + g') + 5Z'$ + 5g') + 

l<p<n-l 



(73) 



> 1, -(qW^I^J = ^n-i((5T') - (5r)|^) + 

5„,i((5f ) - m|^) + ((5T)o,„-i - (5T)o,„-i|^) + 
E Un-p((^')op - (nop|^) +^„-p-i((5T')op - (5r')opi^)). (74) 

l<p<n-l ^ ' ^ ' 

3.5 Perturbation series of a bilinear form of the fluxes 

I seek the power series of (DT)*, where D is any endomorphism on T, hke 
D = d/dz, such that SD = 0. 



Remark 5. By / l5i| , \6^ ), the operator (D.) depends implicitly on T , through 
the fluxes, so that, m general S{{DT), 6 {D.),6T)> A{DT) ^ {6DT). 

I evaluate A{DT) by applying twice the formula for the perturbation of a 
bilinear product (pOj); after rearranging, I obtain the perturbative expression 

A{DT) - {5DT) = 

{mJ*^^ + DQ^*\6^) + {6^^DL*^ + DQ*) + {S^^DL*6^). (75) 



I assume that D is linear. Applying D and (^) on (^31) , noticing that D 
commutes, by linearity, with 6 and ', I obtain that (|6lD, where T — ^ DT, still 
holds, giving the power series of DT *. This power series, and that of {U_,z)* 
(|59|), are inserted into (|75D, yielding the power series 



5i{DT),5{D.),e6T)> A{DT) =J2{DTUT,ST)e^, (76) 

n=l 

the coefficients of which are determined by identification. Eventually, 

oo 

{{DT),6{D.),e6Ty > (DT)* = Y,{DT)n(T, 6T)e\ (77) 



n=0 
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\/n > 0, (DT)„(T, 5T) = zjDT') + z^^^DST') + {DT^, + D5T^_,) + 

Op+pO ~l~ ^n— p-i{DST')op+po + {DTn-p + -D(5T„_p_i)op+po) + 

l<p<n 

^ ^ (^n-pi-p2 (-^-^ )piP2~^^n-pi-p2-l{^^'^ ) PiP2~^ {^'^n-pi-p2~^^^Tn-pj^-p2-l) P1P2) ■ 



l<pi <n— 1 
l<P2<"-Pi 



(78) 



In (|78D, ^„ is known and need not be eliminated, which allows to group 
terms nicely: 



\/n > 0, (DT)„(T, 6T) = 

^ ] (:5n-pi-p2 (-^-^ )piP2"'":^n-pi-p2-l '^-^'^-^ )piP2 + (-DT„_pj_p2 +-0'^^n-pi-p2-l)piP2 ) • 



0<pi,p2<n 
Pi+P2<n 



First and second order results are 

(DT), = z.iDT') + {D6T) + (DT)oi+io, 

(Z}T)2 = z^iDT') + (DT2) +^i((D5T') + (DT')oi+io) + 

(D5r)oi+io + {DT)o2+2o + {DT)n. 

4 Weight definition and observability 

4.1 The balance equation and weight definition 

For an infinitesimal perturbation of the constrained system, (^) becomes: 

= dRridT). (80) 
Applying ( P5| ) to d{T{z)), then the constraint, with (^^1): 

d(T(2)) = T'(2)dz + dT{zy, dT = rdz + dT. (81) 



(79) 



I introduce (|81|) into (|80D , use (|6^) , where D ^ d, to obtain the balance 
equation 

= {T')dz+ {dT). (82) 
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(T') is the 'differential weight' of the exciting variable z. ( p2D is indeed 
equivalent to (|69|) . 

From the chain rule, the differential weight is just the derivative of the 
unconstrained gauge output (at constant T), with respect to z, hence an 
observable: 



(T') = {R o T)'U) . (83) 



I show up the exciting variable Zi, driving 6T, while the original control 
variable is noted Z2- 

3T2,{T + 6T){z2)=T2{zi,Z2). (84) 

For a given exciting perturbation 6T, (the function) T2 remains constant. 

Remark 6. An experiment that is not reproducible is still represented by a 
constant T2, taking for Zi the quasi- stationary time. 

I assume for T2(., 2:2) the same regularity properties as for T2(zi, .). Thus, 
the variables (^1,^2) are exchangeable, and either may serve as the control 
variable. I define an exchange operator: 

ET2{z2,z^)=T2{z,,Z2). (85) 

The constraint operator applies to each partial function T2(., ^2), ^2(-2i, •)■ 
For example, with Z2 as the control variable, 

u{T2{Zi, .)) = u{T2{Zi,Z2)),Z2{T2{zi, .)) t> Z^. 

The functions Z2 ^ z^(T2{., Z2)) and zi 1— > Z2{T2{zi, .)) are inverse to each 
other and the constrained value of u{T2{zi, Z2)) does not depend on the choice 
of the control variable. For all {zi,Z2), such that R{T2{zi, Z2)) = Rq, 

(Zl.Z^) = iZi,Z2) = iz^,Z2), 

u{T2{zi, .)) = u{T2{zi, Z2)) = m(T2(., Z2)). 

With ([8^), the balance equation (|82D (with Z2 as the control variable) 
takes the nearly symmetric form 

= {diT2{.,Z^))dZi + {d2T2{zi, .))dZ2. (86) 

The differential weight of zi (controlled by Z2) is 

wi(T2,zi) = (5iT2(^i,.)), (87) 
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and symmetrically, so that (^) becomes 

= Wi(T2, zi)dzi + W2{T2, Z2)dz2. (88) 
I define the (integral) weight (of Zi), 

Zi(T2,zi) = [ ' wi{T2, z[)dz[. (89) 



The function Zi(T2, .) is the (non-linear) weight scale. By integration of (pSj), 
the sum of weight perturbations, over both variables, is zero: 

Z,{T2, Zi) + Z2{T2,Z2{T2{Z^, .))) - Z2{T2,Z,{T2{0, .))) = 0. (90) 



(|83| , |90|) are the basis of the weighing method, announced by 



4.2 Weight perturbation theory 

From law invariance (1^) (on functionals) : 

= Swi = SZi. (91) 
A general differential weight perturbation is 

Siwi{T2, zi), 0, ((5T2, Szi)) > Awi. (92) 

By definition (|7D, 

Awi = S{{diT2izi,.)),S{di.),S{T2{z,,.))) 
which is just ([76D, where 

T ^ T2{zi, .)-D-^ di;' -> ^2; e 1. (93) 



Using the properties of definite integral in ([89|), an arbitrary weight per- 
turbation is 

<5(Zi(r2, zi), 0, (5(T2, zi))= r wiiT2, zi)dz[ + I ' 5(w7i(T2, z[), 0, ((5T2, 0))dz;. 

(94) 

The last integrand is just (|92D , where 
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( p2D has another important apphcation. I write ( pQ]) as 

Zi{T2,zi) = [ \w,{T2,0),0,{0,z[)ydz[ = wi{T2,0)zi+ [ ' 6{wi{T2,0),0, {0, zi)). 
Jo Jo 

(95) 

The last integrand is just (|92D, where 

zi 0; 6T2 — 0; Szi — > zi. 

Therefore, the weight scale itself is obtained as a perturbation series, which 
is, more precisely, ([76l), where 

6T ^ T2{z[, .) - T2(0, .); T T2(0, .); ^ ^i/ ^ ^2; e ^ 1. 

I complete the determination of the weight scale only in the case of a 
linear exciting variable: dfiT2 = 0, equivalent to 

T2{z,, Z2) = T{z2) + z^5T{z2); 6T{Z2) = d,T2. (96) 

With (0), the remote control condition (^) becomes ^^2^2 = 0. In (|95D , z[ 
is used as the perturbation variable: 



[Wi 



{T2,0),0,{0,z[)r = {{d,T2),6{d,.),z[6Ty>{6Ty = Y,{6T)^zT, 



n=0 



and the coefficients ((5T)„ are obtained from ( [78i where 

e z[; DT ^ 6T; D6T ^ : 

-p)op+po) + 

l<p<n 

^ ] {^n-pi-p2^^'^ ) P1P2 ~^ {^'^n~pi-p2)piP2)l 



l<pi<n-l 



Vn > 0, (5T)„ = {^n-pr~p2 (^^')pi,P2 + (^^n-pi-pa) 



PlP2^ 



(97) 



0<Pl,P2<" 
Pl+P2<n 



First and second order results are 

(5T)i = ^i(5T') + (5T)oi+io, 

(5T)2 = ^2('^^') + (^^02)+a(5T')oi+io + (5T)o2+20+(5T)n. 

Evaluating the integral in (195|) , and with 2:1 — > e, to remember that the 
exciting variable is linear. 



Zi(T2,e) = X](5T), 



n=0 



n + 1 



(9^ 
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4.3 Operator perturbation weighing 



I define an ideal instrument by its parameters T, consisting of linear functions 
of the control variable (there is no interest in non-linear ideal control); the 
measured object is a perturbation 6T, with remote control. The combination 
of the instrument and the measured object is represented by the functional 
parameters T2, 



T2(e, z) = T{z) + e<5T, z(T2(0, .)) = 0. 



(99) 



Observables are not available for T2, because of its ideal part T, but on 
some approximate realization Tg*. 5T2 = Tj* — T2 is the realization error, 
meaning it is impossible to realize exactly the ideal experiment. From (|90D , 



Zi(r2,e) = Z2(T2,z(T2(e,.)) - Z2(T2,0) ^ 

^2(T2*,z(T2*(e,.))-Z2(T2*,0). (100) 



The error analysis is left for the next section. From (p9| , p3D 



Z,iT;,z{T;,e))-Z2iT;,0) 



d2RoT;{e{T;,z),z)dz, (101) 



which I write, more lightly, by impliciting T2*: 



^2(^(e))-^2(0) 



d2R*{t{z),z)dz 



(102) 



As explained in section ^?T| , the functions e and z (at constant T2) are inverse 
to each other. 

explains how to map observables onto weight, by universal opera- 



tions, like summation, independent of the system parameters. ( pSf ) expresses 
what is weight and how to to compute weight, by universal operations, from 
the system parameters, representing the instrument and the measured object. 
The function Zi{T2, .) is a power series, bijectively related to its coefficients 
{6T)n- As control in T2 is linear and remote, ( P?] ) becomes 



Vn>o,(5r)„= Yl (^^) 



(103) 



0<pi,P2<n 
Pl+P2=n 



In particular, from (|62| , pT]) , 

(5T)oo = (5T)o = {6T) = ($t|M) + {^^\SQ) + {6Q^^). 
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There are simplifications also in the recursion relations (0, |6^, |68| ). 

If the sequence of perturbed fluxes is a basis (and adjointly) and 
= 6Q = SQ\ then 6L is completely determined (except for a scalar factor) 
by its matrix elements {5T)p^p,^, defined by (|63D. From ( |103| ), {5T)n is the 
finite sum, invariant by transposition, of the diagonal p {5T)n-p,p- The 
{5T)n do not fully determine (5L, except, for example, if 5L has only one 
coefficient on each diagonal. 

If 7i is not separable (has no basis §3.7, p. 95]), then 5L cannot be 
obtained from observables. We must content ourselves with weight, which 
does not completely determine 5L, but can be used to check postulated 5L, 
given the perturbed fluxes. Or, knowing (5L, perturbed flux calculations can 
be checked, from measured weights. 

4.4 Error analysis 

From (PBD, the error on the measured weight scale is 

z;{t;, .) - z,{T,, .) = {z,{t;, .) - z,{t,, .)) + 6z,{t;, .). (io4) 

The first term in the r. h. s. is the realization error, that may be evaluated, 
just like any perturbation, from weight perturbation formulae ( p4| , |77| , |78|), 
after the variable exchange Z2(T2, .) = Zi{ET2, .) (|85|) . 

The realization error could be removed by taking T2 = T2, thus attaching 
the weight definition to a real (non-ideal) standard ('old-style' metrology). 

The second term in the r. h. s. of ( |104D , (Tg*, .), is a processing error, 
meaning that observables are not processed exactly as demanded by (|101|) , 
for example, because of discretization errors. Processing errors such that 

3f,Z;{T2,.) = Z2{f{T2),.), (105) 

are reducible, by definition ( |105| ), to realization errors. Not all processing 
errors are reducible. 

Errors in the 1. h. s. of (|100|) are mathematically treated just as errors 
in the r. h. s. (except for variable exchange). Although T is ideally defined, 
one may have to compare weight scales resulting from different definitions 
of T. An error on Zi may also occur in computation, e. g. a power series 
truncation. From (0), an error on Zi, which is not reducible in the sense of 
(|105|) , breaks the law of fiux-to-source operator and constraint linearity. 

The weight scale is actually obtained approximately, and only on a finite 
set of e. Something must be said on error propagation, from the weight scale 
points to its power series coefficients. 
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Firstly, I consider only the effect of discreteness of e. 1 assume that 
Zi{T2, e) is exactly known, but only for + 2, > discrete values of e. An 
interpolating polynomial of degree N+1 can be constructed, giving (approx- 
imately) the coefficients {ST)Q<n<N', higher order coefficients are completely 
undetermined: there is a cut-off between known and unknown coefficients, 
at n = A^. If A^ = (only two values of e are realized), then only {ST)o is 
obtained. 

Secondly, I take into account, not only discreteness, but also errors on 
weight values. Errors propagate non-uniformly to the coefficients (5T)„, the 
error on (5T)„ increasing with n, at given A^, and decreasing with A^ at 
given n. Hence the interest of taking large A^, even if only (5T)o is sought. 
A^ > 1 reveals and allows to correct the shielding error on {ST)q obtained 
with A^ = 0. 

5 Conclusion 

I applied perturbation theory to a linear source problem, allowing for a crit- 
ical limit, the consequences of which were non-linear response, ill-posedness, 
a regularizing constraint, and the functional character of parameters. 

The main hypotheses were the constraint linearity and the existence of 
a stable, spectrally separated, harmonic hyperplane. (In spatially extended 
systems, spectral separation usually follows from boundary conditions.) 

The basic tool was linear (or multi-linear) vector algebra. The pertur- 
bation of functional parameters and the constraint application were delicate 
and required specific notations. Progress was achieved by the thorough use 
of symmetries (the last three being familiar in theoretical physics): 

• duality ('adjointness'), 

• commutation relations (between the constraint operator and variational 
symbols) , 

• exchange (of the exciting and control variables), 

• gauge invariance (results are independent of the gauge output unit). 

Perturbation theory was then used to propose a weight definition and 
measurement method. More precisely, the weight defining functional Zi and 
the weighing functional Z2, needed in (^, were constructed, showing weight 
as a secondary observable, except for a realization error, which was analyzed, 
from the same perturbation analysis. 
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The weighing functional is an integro-differential processor (filter) on 
both the unconstrained (open loop) and constrained (closed loop) system 
responses. 

The weight scale Zi{T2, .) is a power series, whose coefficients are (in the 
simplest case) diagonal sums of the flux-to-source perturbation matrix 6L. 
As this relation is generally not invertible, SL cannot be fully determined. 

The present work is not mathematically complete, and does neither treat 
in detail any particular application. Emphasis was put into discovering and 
solving a new problem, operator weighing. 

The mathematical model and the weighing method may actually apply, 
more or less easily, to high-gain feedback linear amplifiers, with any number 
of degrees of freedom: electronic amplifiers, photomultipliers, nuclear fission 
chain reactors 0, 0] and computational models thereof. 

On the mathematical side, progress may be sought in the formulation of 
the abstract problem, the analytic aspects of its solution (the conditions of 
convergence of power series), and carrying out large symbolic calculations. 

Systems are often studied with the approximation of linear response, 
which fails for non- linear systems, but also, as I pointed out, for a linear 
system (especially quasi-critical). Measurement methods based on the linear 
response approximation lack of an essentially non-linear error analysis, and 
hardly have any quantitative interest. 

A Composition of power series 

The power series of the compound function T o z is sought, z has a power 
series at zero and T has a Taylor series at zq = z{0): 



z — 




n=0 



The multinomial formula. 



CO 




Wn>0,{z{e)-Zor = {J2z,en 



n 



p=i 



is multiplied by T^"'^/n\, summed over n, and the terms of same power in e 
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are gathered: 

Toz{e) = 5^(T o z)„(zo)e", (106) 



oo 



n=l 

/Or) 



(To.). = rE..)n^. 

The sum in the r. h. s. is spht into a term corresponding to the sequence 
<ip = Sp^n, and a remainder T„ (Ti = 0): 

Vn>l,(To2;)„ = ZnT' + Tn{Zi,...Zn-l), (107) 

Vn>2,T„(zi,...z„_i) = ^^^''^ n (108) 

p<n,qp>0,Y^pqp=n l<p<n 

To = T, (109) 
Ti = 0, 

T,iz,) = T"zl/2, (110) 
T3(zi,^2) = T"z^Z2 + T"'z!/3\. 
If T is hnear (|4l|), then 

Vn>2,T„ = 0. (Ill) 
Proo/. If T is linear and T^^ ) in ( [10^ ) is not zero, then 

qp = l] 3!m <n,qp = 5p^rn\ ^ pq^ = m = n, 

l<p<n l<p<n 

which is inconsistent. □ 



B Finite-dimensional examples 

The reference gauge output Rq is set to unity. 

I assume, for simplicity, linear and remote control: 

L{z) = A + z'B, 6B = 0, 6L{z') = C = 6A, 
L2(e, z') = L{z') + e6L{z') = A + z'B + eC, 

where A, B, C are constant matrices; A has the eigenvalue zero; sources are 
unexcited (0 = 6Q = SQ'^) and uncontrolled {0 = Q' = Q^'); there is no 
realization error {5T2 = 0). 

I will obtain directly the fluxes, the gauge output, the control variable 



and flux response, and check (pO| , p8| , |102|) . The functional operations, like 
perturbation and constraint, will appear concretely. 
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B.l One dimension 

A = 0, the harmonic subspace is {0}. 

Q 



z'B + eC 



R (e,z) = ' - 



z'B + eC z'B + eC 

With the new variable z = z' + Q'^Q/B, 



Q 



<^*(e,z) 



R*(e,z) 



Q^Q - zB - eC 
Q^Q-zB-tC 



From (13, H, [1021) : 



QtQ 



{Q^Q-zB-eCy ' ' Q^Q' 
Z,{z{e)) - Z,{0) = j^'^^\,R*[eSz),z)dz = ^ = -Z,{e). 

5L/{Q'^Q) is the secondary observable, resulting from the weighing method. 

B.2 Two dimensions, linear control variable response 

The complementary matrix (or transposed comatrix) A of a matrix A is 



/ an 


ai2 




1^21 


^22 , 





—021 Oil 



AA = AA = detiA)l,A = A,A^ = A^. (112) 

The operator A Ais linear and the non-linearity of A i— A~^ = det(A)^^y4 
lays fully in the inverse determinant. 
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Let the 'codeterminant' A* B oi two matrices A, B be 
A*B 



ail 


ai2 


+ 


bn 


bl2 


hi 


b22 


021 


0-22 



Properties: * is a bilinear product, symmetric and, for all A, C, 

A* A = 2det A, A* C = AUc\ 
AC A ={C* A) A - det{A)C. 



(113) 



The reference flux-to-source operator has a stable subspace and the eigen- 
value scale is arbitrary. Therefore, it is no restriction to take 



A 




-1 



Property: for A defined by ( |114|) and for all B, C, 

ACB + BCA ={C* B)A + buC - cuB. 



(114) 



(115) 



Let Q 



and adjointly 



- det L2(e, z') = - det{B)z'^ - det{C)e^ -B* Cz'e + bnz' + cne, 



- det L2$*(e, z') = (A + z'B + eC)Q, 
-detL2R*{e,z') = -q\qi + {Q^BQ)z' + {Q^\CQ)e, 
- det L2 (e) = -qUi + {Q^\BQ)z[ + {Q^\CQ)e. 
I assume, for simplicity, bu ^ and 

= det 5 = det C = 5 * C, 



(116) 
(117) 

(118) 
(119) 



which makes the control variable response linear. ( |119| ) still leaves five degrees 
of freedom in -B, C, out of eight generally. 
With the new variable 



z = z 



«!, ai 



mBQ)-bii 
and (|1T9D, (|116|, |T3, [lT|) yield 

-detL2(e, z) = buZ + cue + huUi, 

{Q^BQ)z + {QMCQ)e + aiiQ^BQ) - q[qi 



R*(e,z) 



zie 



a2e, a2 



biiz + Cue + buUi 

mc_Q)-cii 

mBQ)-bii 



(120) 
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$*(e) = {A + aiB + e{C + a2B))Q, (121) 

1 — 

1 _bn{Q^CQ)-cn{Q'\BQ) 

«3 = 7 >«4 = 1 , 

Oiiq[qi 

aia4 + 02 = -7^- (122) 

Oil 

To avoid the determination of primitives, I will rather check the differen- 
tial form of®: 

-($t*|C$*)(e) = 52i?*(e,z)^(e). (123) 

de 

The r. h. s. of ( |123| ) is evaluated by applying on ( |12C1| ) 

V(A,det A ^ 0),V(^, = 1),- [z) = -— : 

a2i^ + 022 dz a2iz + 022 det A 

dz. , W\BQ)-bu)i{Q^CQ)-c,,) 1 



00 



(124) 

The 1. h. s. of ( |123| ) involves the differential weight of the exciting 
variable: 

00 

($t*|c$*)(e) = w,{e) = Z[{e) = ^(5r„e". 

n=0 

The 6Tn are given by (|103|) , which appears directly by substituting in ($^*|C$*) 
the power series of the perturbed flux ( |121D (and adjointly). The flux and 
weight power series coefficients are 

$0 = a3(A + aiB)Q, $J = 03(1^ + aiB^)Q\ 



$1 = asia^A - ^B + C)Q, = a^^a^A^ - ^B^ + C^)Q^ 
On Oil 



6To = {^l\C%), 

Vn > 1, 5T„ = <-^((<l>S|C$i) + (<fI|C$o)) + a2-\n - l)(<l'I|C<l>i). (125) 
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STo = alQ^ACA + ai{ACB + BCA) + aiBCB)Q. 



From (|113| , |115| , |119| ), and after a few lines of algebra, 



qlqihi 



= aliQ^ia^A - ^B + C)Cia^A - ^B + C)Q) 

Oil Oil 

{Q^\{alACA+ {^)^BCB - a^^iACB + BCA))Q) = 0. 

Oil Oil 

(127) 



Taking ( [I27| ) in ( |125[ ) shows that wi is a geometric series. The ratio 04 and 
the first term given by ( |126| ) are consistent with ( |124] , |123|) . The radius of 



convergence is a^^. Other results can be checked: the spectral properties of 
L{z), the recursion relations on perturbation series coefficients . . . 
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